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Abstract 

A model for the longitudinal structure function Fl in the region of low x and low Q 2 is 
discussed. It is constructed using the hr factorization theorem and a photon-gluon fusion 
mechanism suitably extrapolated to the region of low Q 2 . A phenomenological model 
for higher twist is presented, which is based on the assumption that the contribution of 
quarks having limited transverse momenta is dominated by the soft pomeron exchange 
mechanism. 
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Introduction 



The longitudinal structure function Fl corresponds to the interaction of longitudinally polarized 
photons in the one-photon-exchange approximation of deep inelastic lepton-nucleon scattering. 
The experimental determination of Fl is difficult since it usually requires measurement of the 
energy dependence of this process for fixed x and Q 2 . The structure function Fl vanishes 
within the naive parton model [0] and also in the leading logarithmic approximation of QCD. 
In the next-to-leading order approximation the structure function Fl acquires a leading twist 
contribution proportional to a s (Q 2 ). The main process contributing to Fl at low x is photon- 
gluon fusion, #7 — > qq. The longitudinal structure function Fl is therefore mainly sensitive to 
the gluon distribution. 

In practice it is the ratio R(x,Q 2 ) = Fl(x,Q 2 )/Ft(x,Q 2 ) which is determined experimen- 
tally (where Fj-(x,Q 2 ) is transverse structure function). Current measurements of R(x,Q 2 ) 
(SLAC, NMC at CERN || ) are very scarce, especially at low x. In the region of large Q 2 the 
structure function Fl is fairly well described by perturbative QCD, while at lower values of Q 2 
it is not well understood mainly because of the non-perturbative contributions. A better under- 
standing of Fl and R in this region is important in order to implement the radiative corrections 
needed for the extraction of the Fi{F^) structure functions from deep inelastic scattering mea- 
surements. It is also crucial for the reliable extraction of the spin dependent structure function 
gi(x,Q 2 ) |J from experimental data, since the uncertainties in R are one of the main sources 
of systematic errors. 

In this paper we present a model for the structure function Fl which provides extrapolation 
into the region of low Q 2 . We also propose a model for the higher twist term. It is based on the 
assumption that the contribution of quarks having limited transverse momenta is dominated 
by soft pomeron exchange. We introduce an effective coupling which can be determined from 
the background contribution to the structure function F T ||. Calculations for R(x,Q 2 ) in the 
region of low Q 2 are also presented. This paper is based on the results presented in [||. 



hr factorization vs collinear factorization 

The off-shell approximation is closely connected with the &t (or high energy) factorization 
theorem || [7| which can be presented in the following way, 

Fl(x, Q 2 ) =f d 4i #X 79 ( * , 4, Q 2 )f(*\ k 2 T ). (1) 

The dominant diagram for the Fl is shown on Fig.l, where x', kx are the longitudinal fraction 
of proton's momentum and the transverse momentum of the gluon g respectively. The function 
F2 9 (^7, k 2 -., Q 2 ) can be regarded as the longitudinal structure function of the off-shell gluon with 
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virtuality fcf., and f(x, k^) is the unintegrated gluon distribution which in principle should be 
obtained from the BFKL equation || [|. In formula ([I]) there is double convolution over the 
variables x' and k T , and therefore it is called the off-shell approach. On the other hand the 
well known renormalization group analysis is basically connected with collinear factorization 
formula, 

F(x, Q 2 ) = W 1 ^-Qix/x', a s , Q 2 /fi 2 )x'h t (x', fi 2 ) . (2) 

i 

The hi(x', /i 2 ) are the integrated parton distribution functions which satisfy conventional Altarelli- 
Parisi equations Jl0|| and C{ are the coefficient functions which can be calculated perturbatively. 
The summation is performed over different species of partons and fi 2 is a factorization scale 
which usually is taken to be Q 2 in case of deep inelastic scattering. 

It has been shown in ||, [TTJ] that the leading twist part of the k T factorization formula can 
be recast into the collinear form (fj) where all small x' (x/x') effects are resummed to all orders 
in a s ln(x') (a s ln(x/x')) in the splitting (and coefficient) functions. For instance the moments 
of the coefficient functions Cl{oJ-,ol s ) 

C l (lu, Q 2 , a s (Q 2 )) = jf -z"C L { Zl Q 2 , a s (Q 2 )) (3) 

are related to the double moment function F] 9 ^,^) by 

C L (u,Q 2 ,a s (Q 2 )) = lgg {a s /u)Fl\u )1 = 7 99 (a») (4) 
where a s = N c a s /n, and 

1 r l/2+ioo _ /0 2 \ 7 

F2 9 (u,Q 2 ) = — d 7 F2>, 7 ) ^ 

2m Ji/2-ioD \ k 1 ) 

F2 9 (^Q 2 )= C-x^Fl^x^Q 2 ). (5) 

JO X 

The quantity j gg (a s /uj) is the gluon anomalous dimension related to the appropriate splitting 
function by 

^ gg (a s /u))= dzz^P gg {z,a s ) . (6) 



Knowing the expansion of F] 9 in terms of the anomalous dimension and the expansion of 7 95 
in powers of a s /u, one can obtain the following formula for the coefficient function Cl(z, a s ) 
(in case of F L ), 

C L (z, a.) = a s [Cf\z) + a s C£\z) + a^n(Vf( Z ) + o? s ln 2 (^)Cf {z) + . . .] (7) 

where as z 0, Cf(z) -» C$(z) -> canst ± 0, % > 1. 

If instead of the BFKL solution in (ED), we take 



f(x,k 



2 



dxg(x, Q 2 ) 
dlnQ 2 



2 



(8) 



where xg(x, Q 2 ) is a solution of the Altarelli-Parisi equations without small x corrections in 
P gg but preserving the original phase space (no kx ordering), then only the first three non- 
trivial corrections in small x are included in the expansion for C^. This approximation seems 
to be reliable, but not of course at very low values of x where the BFKL effects could play an 
important role. In our calculations we consider x > 10~ 4 . 



The photon - gluon fusion model for Fl. 

Here we show how after an expanding the integrand in the off-shell formula (|XJ) around k\ = 



one can obtain the on-shell approximation 12] with additional power corrections. If we use 



the kjp factorization theorem to calculate Fl, then we obtain the following form for Fl at low 
x and large Q 2 : 

7T 2 



F L (x,g 2 )=2^^e^(x,g 2 ) (9) 



where 



Q 2 )= l^f £ d/3 J d 2 K' T a s (Q 2 )(3 2 (l - P) 2 \ ^- - ^- ] f(x', 4) (10) 



and 



D lq = 4 + (3(1 - (3)Q 2 + m\ 
D 2q = {kt - k T ) 2 + /3{1 - (3)Q 2 + m\ . (11) 
The transverse momentum k' t and the variable x' are defined as follows: 

4 = k t -(1- P)k T (12) 
/ K' 2 + m 2 q k 2 \ 

X=x \ + w^m + ^) • (13) 

The integration over x' is implicit in (|T0|). One can make it explicit by changing the variables 
from k'j, to x' using the above formula (|13|). The integration limits in ( |10D are additionally 
constrained by the condition: 

x' < 1 . (14) 

The condition x' > x demanded by the integration limit is automatically satisfied, see (|T3D. In 
the on-shell approximation of the formulae d^,|l(3) one expands integrands of the corresponding 
integrals around = and retains the leading term. The integrals I q can then be expressed 
in terms of the conventional integrated gluon distribution g(y, Q 2 ) 

I q (x, Q 2 ) = tt jf 1 dp J d4a s {Q 2 )p 2 {l - P) 2 ^yg(y, Q 2 ) , (15) 
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where now 



and 



y — x I 1 + 



K 



12 



P(l-(3)Q 2 



(16) 



J2 



(17) 

to the low Q 2 region we 
To do this we substitute 



, T + 0(1 - P)Q 2 + m\ . 

In order to extrapolate formula @ (with I q given by ( |T0|) or ([15|) 
have to freeze the evolution of g(y, Q 2 ) and the argument of a s {Q 2 ) 
Q 2 + Am 2 instead of Q 2 as the argument of a s and of yg(y, Q 2 ). 

This substitution may be justified by analyticity arguments i.e. we want I q to have threshold 
singularities for Q 2 < —Am 2 . Moreover, for heavy quarks and for moderately large values of 
Q 2 , it is the heavy mass squared and not Q 2 which should be taken as the scale of a s and 
of the parton distributions. It should be noted that after these modifications the structure 
function Fl can be continued down to Q 2 = 0, respecting the kinematical constraint Fl — > Q 4 
as Q 2 — > 0. The formula ( 15 ) can be rearranged to give the Altarelli-Martinelli [13] integral 
over y but with the modified integration limits and additional power corrections. First we use 
y and /3 as the integration variables instead of k'^ and f3. From (|17D we find 



(3{1-(3)Q 



2,y_ 



D n 



(3{1-(3)Q 



2 y 

x 



From eq. fll§|) we can deduce various integration limits since 



(3(1 - p)Q'(^ - 1) - m >0 
x H 



and since 1 > /3 > 0, 



j 



From inequality (|21[) we obtain the lower limit for integration over y 

Am 2 

y>x(l+ g ^ 



Q 2 



and of course, 



(3(1 - (3) > 



Q 2 (I-i) 

It is convenient to make another change of the integration variables: 

1 



P 



+ X 



The inequality ([H]) gives the following limits for the variable A 



1 m z q 
Q 2 (y--1) 



< A < 



mi 



\ 4 Q 2 (l - 1) 



(18) 
(19) 

(20) 
(21) 

(22) 
(23) 

(24) 
(25) 
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Finally we get the following representation for F^: 



where 



ft,- 



dy ( x 



n •/*« y \y, 



1 - 



m = o^f'dyfA ln 

7T Js q y \yj 



yj\ Q 2 (y-x 



4m « X y 9 (y,Q>) 



l + A 1 



4m^x 



Am^x 

r(y-x) ) 



yg(y,Q 2 ) 



and 



x q = xll + -^- 



It should be noted that expressions (p6| . |27|j28|) reduce to the following formula 



F L (x,Q' 



ot s {Q 2 ) r 1 dyfx s 
I* y \y. 



71 



X 



l--jyg(y,Q 2 



(26) 



(27) 
(28) 

(29) 



(30) 



given by in the Q 2 — > oo limit. The expression -j^-J^ represents the power correction 
enhanced by the logarithmic factor InQ 2 /m 2 q . The modification of the lower limit becomes 
obvious, if we consider the energy-momentum conservation constraint for the subprocess 75 — > 
qq. 



A phenomenological model for the higher— twist contribu- 
tion to Fi 

The kinematical constraint restricts the behaviour of Fl in the region of small values of Q 2 . 
Analysing the longitudinal cross section for deep inelastic scattering one concludes that in order 
to remove the possible singularities of this observable at Q 2 — ► 0, (i.e. o L ~ F L /Q 2 ), F L should 
behave as Q 4 . This implies that ~ Q 2 , so it vanishes in the real photoproduction limit as 
expected. Phenomenological analysis of the quantity R in the region of moderately large Q 2 
( Q 2 > lGeV 2 ) implies that the higher twist effects (i.e. terms that behave like l/Q 2 as 
Q 2 — > 00) should play an important role ]14] . 



In this paragraph we introduce a phenomenological model for the higher twist term which 
satisfies the constraint Fl ~ Q 4 as Q 2 — > 0. The idea is to treat the contributions from low and 
high quark momenta in a different way. We divide the integration over (in (|T0|) for the off- 



shell model (or in (27 ) , (28 ) for the on-shell model ) into two parts - the region < k't < k'q T , 
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and k't > k'q T ~ where k' 2 t is an arbitrary cut-off, chosen to be of order lGeV 2 . The region of 



high momenta is treated in the usual way (using off-shell (|9|,|10D or on-shell (|2^|2^J2^) formula) 
whereas for low momenta we assume the existence of an effective coupling constant: 



a s (Q 2 )zg(z,Q 2 



A . 



(31) 



The constant A is not a free parameter. We can estimate it from F 2 assuming that the non- 
perturbative contribution to F 2 in the scaling region also comes from the region of low values of 
k'j- and is controlled by the same parameter. The contribution to Fl from the region < Kt < 
k'q T is calculated using on-shell formula, making the substitution (|31|) and integrating over k®, 



P)Q 2 + 



I (3(1 - (3)Q 2 + m 2 q + 3k' 2 t ' 
'6 [P(l-P)Q 2 + m 2 + K> 2 T }\ 



(32) 



This formula vanishes like 1/Q 2 for large values of Q 2 , therefore it has the form of a higher 
twist term. It also obeys the kinematical constraint for the structure function Fl , i.e. vanishes 
like Q 4 for Q 2 -> 0. 

To evaluate the constant A we consider the corresponding formula for Ft, 



F t (x,Q 2 



2E' 



2 tf 

Hti 2 Jo 



dp d k!toi s (k!t + "\>) x 



x Up 2 + (1 - P) 2 } 



U \q 



D lq D 2g 



+ 



f(x',k, 



(33) 



where k' t , D lg and D 2q are defined by eqs. (|TT| , [12]). If we rewrite the integrands to be 
symmetric in variables and k t — k T we obtain, 



u \q 



K T ■ (n T — 

D lq D 2q 

2 

mi 



D 

m 



lq 

2 r 



D 



2q 



D( q D lq D 2q 2 



-1 2 



l D lq 



D 



2qj 



(34) 



If we now expand the right hand side of the eq . (|34]) in powers of k T /Q 2 ( on-shell approximation) 
we obtain the following formula, 



k,t — k'i 



1 2 



D 



lq 



D 



2 q 



kx 2k't{k't • kq 



iDq 



D 2 

v 9 



4cos 2 (f)K' T ] 4cos 2 4>k'^ 



D 3 
q 



F> A 

q 



(35) 
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where we have dropped terms containing odd powers of coscf) which give vanishing contribution 
after performing intergation over d(j). Then Ft takes the form 



F T (x,Q 2 ) = 2A£^ I' d(5 f <T dn% 
q H 4tt Jo Jo 



x {_^ + (1 _^__^ + _^ + _^j. (36 ) 

Performing the integration over k't, substituting (|3~I~D and taking the limit Q 2 — > oo , we get 
the following form for the formula for the structure function Fl, 



a - 2nF - (37) 

{2 m q 1 K 0T 1 

The value for F T 9 is taken from the phenomenological model for Ft |§] where Fy is assumed to 
have the form, 

F 2 ~ F T = F T 9 + F^ FKL (38) 

where F^. 9 is the background originating from non-perturbative region and its value is fixed by 
the data to be around 0.4. Assuming k' 2 t = lGeV 2 one obtains A ~ 1.96. 



Numerical results 

In this section we present the numerical analysis for the structure function F^. We shall also 
determine the ratio R and compare it with SLAC parametrisation |15|]. We include u,d,s,c 
quark contributions, with quark masses 0.35, 0.35, 0.5, 1.5 GeV respectively. We have used both 
GRV [0 and MRS (A) |L7) parton distributions. In GRV the LO gluons and coupling constant 



were used. The LO coupling constant has the following form, 

a s (Q 2 ) 2 



2vr ' p l n 



(39) 



A- 



where (3q = 11 — 2JV//3, Nf is the number of flavours, and the value of the constant A depends 
on number of active flavours 

Ag^ 4) = 0.232, 0.200 GeV . (40) 

We include the charm quark production threshold i.e. if Q 2 > Am 2 then Nf = 4, otherwise 
A / = 3. 

In order to extrapolate the structure function Fl to low Q 2 we use the method of freezing the 
argument of the coupling and of the gluon distributions as explained before. For the off-shell 



formula we have used the prescription proposed in [18]. In order to regularise the integration 



over k T we have introduced the cut-off h,Q T . The contribution of integration from to /cq T 
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is calculated using the approximate formula (^,^,^) with the substitution yg{y,k1 T ). The 
magnitude of the structure function Fl turns out to depend rather weakly on the chosen cut-off 
k^rp. The dominant contribution to F L is given by the photon-gluon fusion subprocess — > qq 
(so the gluon distribution is the most important at low x) but the contributions from quarks 
have also been included with the modification x — > x(l + Am 2 /Q 2 ) for the lower limit of the 
integration over y. 

AF L (x ; g 2 )=Ee f s(g2 + 4mf) ^ / 1 ^(-) 2 fe(y,Q 2 + 4m l 2 ) + g^Q 2 + 4m l 2 )] (41) 
i 7r 3 Jx q y y 

where x q is given by ( ^9|) . In Fig. 2. we compare our results for massive and massless quarks. 
As can be observed the quark masses play an important role in the region of low Q 2 . But 
even at moderate values of Q 2 the difference is visible. This is because the on-shell model 
(|26| , |27| , |28|) with masses contains the power correction term multiplied by the logarithmic factor, 
~ m 2 / Q 2 ln(Q 2 / m 2 ) . Note also that the difference between the on-shell and off-shell calcula- 
tions is very small. 

In Fig. 3 we present Fl as a function of x, with and without higher twist corrections. We 
have chosen the cut-off parameter k'q T to be l.QGeV 2 , but there is no significant dependence 
on this quantity (within the range of 0.8 < Kq T < 2.0GeV 2 ). Higher twist gives higher values 
for Fl, especially in the region where x is not very small (x > 0.001). On the other hand it 
gives lower values for small x than the model without higher twist. This is due to the fact that 
the contribution coming from low values of the quark momenta is now described by the soft 
pomeron exchange. The difference between these calculations decreases as Q 2 increases. 

In Fig. 4 we present the results for the ratio R(x,Q 2 ) = F L {x,Q 2 )/(F 2 (x,Q 2 ) — F L (x,Q 2 )) 



together with SLAC |15j parametrisation. For both parton sets (GRV and MRS (A)) the struc- 
ture function F 2 has been calculated from the VMD model We also observe the same 
effect as for the longitudinal structure function: the model with higher twist predicts bigger 
values for the ratio R at high x. 

For intermediate values of Q 2 this model coincides with SLAC parametrisation, whereas 
for lower values of Q 2 it vanishes quickly as Q 2 — > 0. It should be noted that the difference 
between the parametrisations in the low Q 2 region are not very large, they both decrease very 
fast with decreasing Q 2 . This suggests that the ratio R and the structure function F L are not 
very sensitive to the differences in parton parametrisations in this region. Because the SLAC 
model is not applicable in the low Q 2 region, using it as an estimate of the ratio R may cause 
significant errors. 



Summary 

We have analysed the proton structure function Fl in the low Q 2 and low x region. We have 
studied the off-shell formula which originates from the high energy factorization theorem and 



8 



we have proved that in the limit when k\ — > it reduces to the on-shell approximation. The 
numerical differences between these two approaches are not very significant. We have also 
studied the structure function Fl in the very low Q 2 region, where we have proposed the 
method of continuing the formulae for the structure functions down to Q 2 = 0. The non-zero 
values of quark masses have strong consequences in the structure function behaviour. Both 
the off-shell and on-shell formulae contain power corrections connected with quark masses (i.e. 
terms like m 2 / Q 2 ln(Q 2 / m 2 ) ). We have also proposed a model for the higher twist term which 
is based on the assumption that the contribution of quarks having limited transverse momenta 
is dominated by soft pomeron exchange. This results in change of the x dependence of the 
structure function F L . We have also performed calculations for R. The model for the ratio R 
enables us to estimate this quantity in the low Q 2 region. It is particularly important because 
the uncertainty in the ratio R is the main source of systematic errors when extracting other 
structure functions from experimental data. 
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Figure captions 

Fig.l The diagramatic representation of the fcr factorization formula for Fl. 

Fig. 2 F L plotted as a function of x for different values of Q 2 . We use GRV partons [pi. Solid 



line: off-shell approximation with k\ T = 0.23GeV 2 , dashed line: off-shell approximation with 
k\ T = l.OGeV 2 , dotted line: on-shell approximation, dashed-dotted line on-shell approxima- 
tion with massless quarks u,d,s. 

Fig. 3 Fi as a function of x for different values of Q 2 calculated in the on-shell approximation: 
comparison of higher twist effect, partons distributions are obtained using the GRV parametri- 
sation. Solid and dotted lines correspond to the results without and with the higher twist term 
respectively. 

Fig. 4 The ratio R plotted as a function of Q 2 for different values of x with the higher twist effect 
included and compared with the SLAC parametrisation (thick solid lines). Solid and dashed 
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lines correspond to the off-shell and on-shell approximation with the GRV parametrisation of 
the parton distributions |]I6| . The dotted and dashed-dotted lines correspond to off-shell and 
on-shell approximation calculated using MRS (A) partons [ITT . 
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Fig.l 



Fig. 2: F L (x,Q 2 ) 
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Fig. 3 F L (x,Q 2 ) 




Fig. 4: R(x,Q 2 ) 




